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Abstract 

^ ' Dynamical aspects of quantum Brownian motion in a low temper- 

r , ature environment are investigated. We give a systematic calculation 

of quantum entanglement among two Brownian oscillators without in- 
voking Born-Markov approximation widely used for the study of open 
systems. Our approach is suitable to probe short time dynamics at 
cold temperatures where many experiments on quantum information 
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1 Introduction 



Chemical systems such as hquid NMR have been actively used for quantum 
information processing experiments[l]. Hybrid devices based on molecular 
ensembles on solid state circuits provide many advantages as quantum pro- 
cessors such as long lasting memory, good controllability and scalability [2]. 

Implementing quantum information processing devices requires precise 
knowledge of quantum open systems. The detailed study of quantum open 
systems are often prevented by a limitation of analytical or numerical re- 
sources to probe systems with large degrees of freedom. The first principle 
calculation of open quantum systems based on projection operators[3, 4] 
or influence functionals[5] has been developed and applied to numerous 
problems[6, 7]. Reduced dynamics obtained by these methods follows non- 
Markovian evolution carrying the memory of environment. Solving the equa- 
tions of motion obtained by these method directly is generally difficult to 
deal with as they are given by integro-differential equations. 

Various approximation schemes have been used to simplify the exact evo- 
lution to study open system evolution. Most commonly used approximation 
is Born approximation in which the system-bath interaction is treated per- 
turbatively. However, naive Born approximation applied to quantum open 
systems does not guarantee the positivity of the density matrix and requires 
a caution particularly when evaluating quantities sensitive to this aspect 
such as quantum entanglement. 

For a weak system-bath interaction, bath is disturbed by the interaction 
only in a small amount and tends to reset to the original state in a short 
time compared to the system time scale. The correlation among the bath 
variables may become negligible compared to the system time scale. In a 
Born-Markov approximation, the bath correlation is ignored. At sufficiently 
low temperature, however, the bath memory can persist for a long time 
and this approximation is expected to break down. Regardless of the tem- 
perature, the noise with long range correlation such as 1/f noise cannot be 
treated under this approximation. 1/f noise is ubiquitous causing a prob- 
lem in implementing a quantum circuit using electronic instruments in solid 
state, ion trap, or hybrid molecular devices. Decoherence in a nonequilib- 
rium bath[8] also shows a peculiar behavior which is significantly different 
from the one under Markov approximation. 

Further approximation such as rotating-wave approximation is commonly 
used with Born-Markov approximation. There we assume weak system-bath 
interaction such that a coupling becomes only among near resonant modes. 
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We further ignore counter rotating terms since they are expected to be av- 
eraged out during a sufficiently long observation time. The rotating-wave 
approximation is thus not suitable for describing short time dynamics in a 
strong coupling regime. 

Quantum entanglement, manifestation of intrinsic nonlocality in quan- 
tum mechanics [9], now became one of the most active research topics in 
quantum information science. It also has been a long historical issue as the 
discrepancy between quantum and classical mechanics is a serious obstacle 
to understand macroscopic classical and quantum mechanics from the view 
of microscopic quantum dynamics[10, 11]. 

Conversely but for the same reason, to manipulate quantum systems 
in order to be useful in the macroscopic world is a formidable task. Recent 
progress in this direction is motivated by rapid progress in quantum technol- 
ogy, where there is a demand for a precise control of small quantum devices. 
Quantum entanglement is considered to be a valuable resource for practical 
applications such as quantum computation and communication[12]. 

Many works have been devoted to clarify the rigorous criteria for entan- 
glement. For continuous variables, the necessary and sufRcicnt criteria for 
entanglement can be given in terms of Peres-Horodecki criteria[13, 14, 15], 
negativity [16], and entanglement of formation [17]. Most of these works are 
devoted to study the static properties of general pure and mixed states. In 
light of information processing, it is also desirable to study the dynamical 
properties based on the realistic models in an open system setting. 

Quantized harmonic oscillators have been playing an important role in 
a history of quantum mechanics. Realization of quantum protocols based 
on continuous variables such as quantum teleportation[18] and quantum key 
distribution[19] show that they also play the similar crucial role in quantum 
information science. Realization of the similar protocols based on solid state 
devices is a highly formidable task since the disturbance from environment 
in solid state systems is much stronger than in optical systems. 

So far the studies of continuous variables for open systems are based 
on the master equation under Born-Markov rotating wave approximation 
such as Lindblad equations [20] or the phenomenological model with high 
temperature Markovian bath[21]. Applications of these methods are limited 
to the high temperature regime or slow dynamics with Ohmic noise, while 
actual solid state implementation operates at low temperatures manipulated 
with fast pulses and often suffers from 1/f noise. 

Most of the solid state implementation of quantum information devices 
operate in a ultra-cold temperature in order to maintain quantum coherence. 
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For instance, typical superconducting qubits operate less than 100[mii'][22]. 
At ultra-cold temperatures, the memory effect in the environment variables 
is no longer negligible. The system-environment interaction in solid state 
devices is typically much stronger than those qubits based on atom-optical 
devices. These conditions are opposite to those where the Born-Markov 
approximation is applicable. Thus we need to develop the method to probe 
the regime beyond the conventional Born-Markov approximation. 

Born and rotating wave approximations are limited to weak coupling 
regime and do not guarantee the positivity of the density matrix evolution. 
Since separability criteria is based on the positivity of the density matrix, a 
caution is required to study entanglement under these approximations. 

In the present paper, we study the precise mechanism of open system en- 
tanglement dynamics based on a Brownian oscillator model[23], where both 
system and environment consist of harmonic oscillators. Since our model is 
formally exactly solvable, it makes us possible to probe the precise dynamics 
of entanglement without conventional approximations. Previous works for 
the study of a two-level system revealed that non-Markovian dynamics that 
plays the major role at low temperatures is highly nontrivial[6, 24]. 

In section 2, we develop systematic tools for our analysis. The reduced 
density matrix in the Wigner representation is calculated for two quantum 
Brownian oscillator model. In section 3, we first make sure that uncertainty 
relations are always satisfied in our model. They are time and temperature 
dependent as our system stays mostly far from equilibrium. Next we make 
partial transpose operation to our density matrix. In the Wigner repre- 
sentation, the partial transpose operation corresponds to the partial mirror 
reflection. We again study uncertainty relations after partial transpose. Ac- 
cording to the Peres-Horodecki-Simon criteria, the violation of uncertainty 
relations after partial mirror reflection can be used as a signature of quantum 
entanglement. Entanglement measures, the negativity and the logarithmic 
negativity, are also calculated for comparison. 

2 General formulation 

We consider a system composed of two harmonic oscillators. Our Lagrangian 
is given by 

^5 = E^-^(^l'^2), (1) 
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where the potential V is assumed to be harmonic: 



2" 



F(i?i,i?2) = E^4 (2) 



Each oscillator variable Rj located at the origin x = couples linearly with 
a scalar field (p via minimal coupling as 

2 

Li = -Y^Xj Rjit)(l)iO,t). (3) 
3=1 

The field Lagrangian is 

L4> = IJdx [(5^0)2 - 771^02] . (4) 

The scalar field (p propagating in one-dimensional space allows a mode de- 
composition: 

<j>{x, t) = j ^^^yf^^^ {6,e--'=*+^^- + 6ie-'=*-^-} . (5) 

In this paper, we study the massless field, then uj^ = 1^1- We consider the 
field ^ as an environment and trace out to obtain dissipative dynamics for 
oscillator variables Rj. 

The Heisenberg equations that Rj satisfy have the form of damped har- 
monic oscillators [7]: 

M,^^^ + M,^)R,{t) -j: f^^ dai,,it,s)dR,{s) ^^^ 

(JjO ^ ^ t/0 CZt' Uo 

where 

otijiit, t') = -2\jXi sin [uJk{t - t')] /ojk (7) 
k 

is an imaginary part of the response function[5] defined as 
aji{t, t') = 2\j\i J2k e-''^'=(*-*') /wk with T,k = I dk/{2Tr). Note that aiji{t, t') 
is antisymmetric in indices: ajji(t,t') = —aiij{t',t). Without any approxi- 
mations, Eq.(6) has a nonlocal form with kernels given by ajji{t,t'). Thus 
the value of Rj at each moment depends on their entire history of the past. 
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In one dimensional space when a ultraviolet cutoff of the field modes 
is brought to infinity, Eq.(6) will be reduced to a local form. Wc write 
71 = Af/Mi,72 = Al/M2,7i2 = AiA2/M2,72i = A1A2/M1. In this case, 
Heisenberg equations of motion will be 



M2R2{t) + M2^lR2{t) + 72M2i?2(i) + 72iMii?i(t) = 0. 



(8) 



We write a pair of solutions of (8) with initial conditions -Ri(O) = -R2(0) = 
and i?i(0) = 1, i?2(0) = as hi{t) and h2,{t). For identical two oscillators 
(Ml = M2 = 1, 0,1 = 0,2 = O) coupled to with equal strength (Ai = A2) 
(hereafter we will drop indices from these parameters), the solutions are 
given by hi{t) = {gi{t) + go{t))/2 and hsit) = (51 (t) - 5o(t))/2, where 



51 (i) 



sin(riri) --Yt J sin(rii) 
'-e ^ and c/o(i) = 



O 



(9) 



are solutions corresponding to two normal modes of a coupled oscillator and 

General solutions with arbitrary initial conditions Rj^ and P^o of coupled 
Heisenberg equations (8) for j = 1, 2 are 

^3 (*) = C'ilj Ri -^10 + C'iij Pi -Plo + (^^3^2 -^20 + C'Pj P2 -f*20 



+ 



+ 



A / dsgx{t - s)^{s), 
Jo 

Cp^R^Rio + Cp.p^PiQ + Cp.ji2R2o + C'p^.p2-f2o 

A / dsg2{t - s)(j){s), 
Jo 



(10) 



where g2 = gi- The expectation value of phase space variables can be 
expressed in a matrix form: 

/ (Ri) \ 
(Pi) 
{R2) 

V (^2) J 

A time evolution matrix C for our solutions in (9) is given by 



C 



1 Rio \ 






CriPi 


CR1R2 


CR1P2 


\ 




^ Rio 


\ 


Pio 




CpiRi 


Cp^p^ 


CpiR2 


CPlPz 






Pio 




R20 




CR2RI 


CR2P1 


CR2R2 


CR2P2 






R20 




\ P20 ) 




\ Cp^Ri 


CP2P1 


CP2R2 


CP2P2 


) 




\ P20 


J 



(11) 



c = 



( fi hi h \ 

f2 h2 fi /14 

h hs fi hi 

V /4 hi f2 /l2 J 



(12) 
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where f2j-i = h2j - h2j{0)gi and f2j = f2j-i for j = 1,2. 

It is convenient for our purpose to introduce the Wigner distribution 
function as 



W{Ri,Pi,R2,P2) = 



(27r)2 J " -"^"^ 2 '"^ 2 '"^ ■ 2 '"^ ' 2 



The characteristic function[27] for the Wigner distribution is given by 



xw{y,t) = IV 



p{0)e 



(13) 



where we defined X = {X1...X4) with ^2j-i = V^Rj, X2j = Pj/V^ for 
j = 1,2 and y = (Y1...I4). We will fix the normalization for each component 
that appeared in C accordingly. The symmetric correlations can be obtained 
from xw{y,t) as 



{{Xi,Xj}) = - 



dYidYj 



\y=o, 



(14) 



where {A,B} = (AB + BA)/2 is an anticommutator. We trace out the 
field (f) in order to obtain the reduced dynamics of the system. With a 
factorized initial condition: p(0) = P5(0) ® p,f){0), xw is also factorized to 
two components as Xwiy,t) = Xwiy ^'^)Xwiy^'^)- '^^^ case, the system 



part is Xwiy^t) = Trs 



, where Xcj{t) are solutions 



of Heisenberg equations with (p = 0. The field characteristic function in (13) 
can also be evaluated exactly. We assume that environment is initially in a 
thermal state with an inverse temperature /3 = 1/T. Its density matrix is 
given as p^{0) = J2k^~^'^'' I ^)(^ I- We obtain 



Xw{y,t) = exp 



exp 



It 



--(F,...n)^ 



/ En ... E 



14 



E44 



\Y, 



■ (15) 



Here 



^2 j-oo rt rt 

T,u(t) = 7^/ (koLoe~'^'^coth.{Pio/2) ds ds'g-At - s)cosu;(s - s')gr(t 
27r ^0 Jo Jo 



where j = (3 + {—iy)/2, are time-dependent (nonequilibrium) fluctuations 
of the system variables, a part induced from environment. Here we intro- 
duced the cutoff frequency A for the field modes. Note that off-diagonal 
correlations "Eji for j I are nonvanishing in general, that is, an interaction 
with a common environment induces an effective interaction between two 
oscillator variables and thus correlation and entanglement between them. 

For an initial Gaussian state with vanishing mean positions and momen- 
tums, (Af(0)) = 0, the system characteristic function also takes a Gaussian 
form: 



Xwiy^t) = exp 



= exp 



(17) 



({XicXic}) ... {{Xci.Xca}) W li Y 
{{XcA,Xci]) ... {{Xca,Xca]) I\Y^ j 



where Xc = {Xci---Xca) satisfy the equations of motion (8) for damped 
harmonic oscillators. (AAf)^(i) are essentially the initial fluctuations of 
system variables shifted by damped oscillatory motion of a coupled harmonic 
oscillator. Combining with the characteristic function for the field, we obtain 



W{X,t) 



(27r)2 (det(A-Y)2(i))V2 
where {l^Xf = {AX)l + E. 



exp 



-\x^{{AXf{t))-'X 



3 Entanglement dynamics of quantum Brownian 
oscillators 

Let us consider a two-mode squeezed state with a squeezing parameter r as 
an initial state [26]. Its correlation matrix is 

(A^)j(o).({^.(o,.^s(o),) = \{:^^^ ITS")' 

In the Wigner representation, the same state can be expressed as 
W{Ri,R2,Pi,P2) = Ae-e'1^^(«i--R2)'+{-Pi+-P2)Vf^]-e-2'-[n(iii+iZ2)2+(Pi-P2)Vf2].(i8) 
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This state can be obtained by acting a squeezing operator e^'^iRi^^-PiRi) 
the vacuum. 

A criteria for separabiUty of a bipartite two-level-system was studied in 
[13]. The necessary and sufficient condition for separability of the density 
matrix is to have only non-negative eigenvalues after partial transpose of one 
of its subsystem. The same criteria is not always sufficient for the bipartite 
system with more than two levels [14]. The extension of this criteria to con- 
tinuous Gaussian variables was first given in [15]. For Gaussian variables, the 
partial transpose of a density matrix in a coordinate representation for one 
oscillator component is equivalent to a mirror reflection of that component in 
the Wigner distribution function. The necessary and sufficient condition for 
a continuous variable quantum state to be separable is that the partial mirror 
reflected state is still a physical quantum state that satisfies the uncertainty 
principle. In the phase space representation, the partial mirror reflection on 
the second variable can be stated as {Ri, Pi, R2, P2) {Ri, Pi, R2, —P2)- 
In terms of X, this can be expressed as the matrix operation by the matrix 
rf = diag(l, 1, 1, —1) as X ^ -qX. It follows that the partial mirror reflection 
transforms the covariance matrix as 

{AXf ^ TiiAXfri^. (19) 

Before we apply the above criteria to our covariance matrix (AAf)^, let us 
make three local Bogoliubov transformations to simplify the form of the co- 
variance matrix. These local transformations do not change the separability 
of the system. First we consider the local orthogonal transformation 

M4=(^^' j eO(2,i?)(g)0(2,i?)cO(4,i?), (20) 

where M2 G 0{2,R) is an orthogonal matrix. Under M4, a symmetric 
matrix of the form 

{AX)^^{{X,X'^})^^^^ (21) 



transforms to 



M,{AX) M, = I ^^^tmT m2DM^ ) ■ 
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By a suitable choice of M4, we can diagonalizc D. Next we make a local 
symplectic transformation with another matrix S/^ G Sp(2, R) Sp(2, R) 
that has a form: 




(23) 



where 5*2 G Sp{2, R) is a symplectic matrix. For a suitable choice of 
<S'2, we can make 52 M2-DMJ S2 to be diagonal with an equal compo- 
nent d. Furthermore another transformation with an orthogonal matrix 
O4 G 0(2, R) 0(2, i?) that has a form: 



04 = 



Oi 
O2 



(24) 



can make (AAf)^ into the following canonical form: 



/ d 


a 
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\ 
b 





(25) 



From the Williamson's theorem [28], there exists a symplectic transfor- 
mation that diagonalizes a positive-definite 4x4 symmetric matrix into the 
following form: 



f Ci \ 

Ci 

C2 

V C2 y 



Although such a symplectic transformation does not preserve the eigenvalue 

spectrum in general, the diagonal components Q for I = 1,2 can be calcu- 
lated. Writing a commutation relation in a matrix form as = iT 
with 



( 



r 



V 








1 



-1 



\ 

1 

0/ 



we construct a real symmetric matrix AA'r(AA')^r^AA'. This matrix has 
an eigenvalue spectrum Q {I = 1,2) [29]. The uncertainty relation can be 
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generalized to a symplectic invariant form (AA:")^ + iT /2 > 0. By changing 
to the diagonahzed form (AA"/))^ above, the uncertain relation is equivalent 
to saying that C,i ^ 1/2 for all I. Although this uncertain relation is invariant 
in arbitrary symplectic transformations, they can change the entanglement 
property. Thus we restrict our transformation to local symplectic transfor- 
mations and use the canonical form of (AAc*)^ in (25) for our analysis. Our 
separability criteria is invariant under local transformations. 

The uncertainty relation expressed by the components in (AAfc)^ are 
given by 

{d + a){d + h)>\ 

{d-a){d-h)>\ ■ ^'""^ 

These are generalizations of the familiar uncertain relations for pure state 
two oscillators to general mixed states. They can be expressed as A^f APf > 
I and AR2AP2 > I in the coordinates that diagonalize the correlation ma- 
trix (AA')2. 




Figure 1: The temporal evolution of the uncertainty function ( in Eqn. 
(28) ) is plotted. The initial condition is a two mode squeezed state with r = 
0.05. Other parameters are 7 = 0.01[GHz], A = 50[GHz], n = 1.0[GHz]. 
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Under the partial transpose in (19), b —b. Thus the conditions for 
separabihty can be written as 

{d + a){d-b)>\ 
{d-a){d + b)>l. 

One can easily see that the two-mode squeezed state wc introduced in (18) 
satisfies the uncertainty relation in Eq.(26). The separability condition (27) 
implies e*^*" > 1 and e"^'' > 1, which only holds if r = 0. Thus this state is 
entangled as long as r 7^ 0. In r — 00 limit, the state becomes the original 
EPR state discussed in [9]. 

In order to make invariance under the local Bogoliubov transformation 
manifest, one can write uncertainty relations (26) explicitly by using the 
symplectic invariants constructed from the determinants of covariances \A\, 
\D\, \{AX)^\ as 

= \D\ + \A\ ± ^{\D\ + \A\r-\{AX)-\ > \. (28) 

In Fig. 1, the temporal behavior of the uncertainty function C,- is plot- 
ted. The initial state is a pure two mode squeezed state and satisfies the 
minimum uncertainty 1/4. As the state becomes mixed, the uncertainty 
increases monotonically in time even for a zero temperature case. At higher 
temperature, the rate of increase is faster. 

Similarly the separability conditions (27) are 

\l = \D\ - \A\ ± ^{\D\-\A\f-\[AXr\ > \. (29) 

Note that the inequalities for C,j^ and A+ in (28) and (29) follow automatically 
from those for ^_ and A_. Thus A_ carries the essential information on 
the separability of quantum states. In Fig. 2, time evolution of the A_ 
is plotted. For an initial coherent state, the uncertainty relation for the 
partial transposed state is always satisfied throughout the whole evolution 
indicating that there is no entanglement. For an initial squeezed state (r = 
0.1), the uncertainty relation is violated initially but eventually satisfied 
indicating that there is a crossover from an entangled to a separable state. 
The asymptotic value of separability seen in A_ appears to be independent 
of the degree of initial squeezing. 

The negativity A/'[16, 30] for quantum Brownian particles can be defined 

as 

N = Mlzl^ (30) 
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Figure 2: The temporal evolution of the uncertainty function ( A_ in Eqn. 
(29) ) after partial transpose is plotted. The initial condition is a two mode 

squeezed state with r = 0.1 for the thick solid curve and a coherent state 
(r = 0) for the thin solid curve. Other parameters are 7 = 0.1[GHz], A = 
50[GHz], n = 1.0[GHz], T = 0. 



where pj is the reduced density matrix after partial transpose. J\f is equal 
to the sum of all negative eigenvalues of and measures how much pj fails 
to be positive. From the Peres criteria, it can be used as a measure of entan- 
glement. It also has a nice property as an entanglement monotone such that 
it does not increase under local operations and classical communications. 
The logarithmic negativity defined as 



Em = logsllp; 



(31) 



also has the similar property. Since diagonalization of (AA')^ brings the 
original state into the thermal state, the partial transposed density matrix 
after the same transformation also has the thermal form that can be written 
as a function of the symplectic invariants X± as 



n 

± 



pA± - 1 
2A± + i;^^l2A± + l 



E 



|n±)(n±| 



(32) 
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For separable states, \± > 1/2. Thus ||p^|| = 1 and Af = Ej^ = follows. 
For entangled states, A_ < 1/2 but > 1/2. The latter follows because 
A+ > C- for 1^1 < (if 1^1 > 0, the state is separable from (26) to (29)). 
WpJW = 1/2A_ follows. Thus both J\f and Ej^j- can be expressed in terms of 
A_ as 



max 



0, 



1 -2A_ 



4A_ 

Em = max[0,-log(2A_)]. 



(33) 



In Fig. 3, the negativity J\f and the logarithmic negativity Ej^ are shown as 




Figure 3: The negativity J\f and the logarithmic negativity Ej\f are shown. 
The initial state is a two-mode squeezed state with r = 0.1. M and Ej^ both 
remain zero in case of the two mode coherent state initial condition r = 0. 
Other parameters are 7 = 0.1[GHz], A = 50[GHz], n = 1.0[GHz], T = 0. 

a function of time. The initial pure two mode squeezed state has the max- 
imum entanglement that decays monotonically in time. They both vanish 
at toE = 0.45[ns] and remain zero. This disentanglement time scale toE 
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measured this way is the same as the time when the uncertainty relation 
for the partial transposed state recovers (the lower curve in Fig. 2) as we 
expect. 

4 Summary 

In this work, we used two quantum Brownian oscillator model to study the 
dynamical aspect of quantum entanglement without Born-Markov approx- 
imation. We studied several different criteria for quantum entanglement. 
The uncertainty function for the two-mode squeezed state under partial 
transpose initially violates the uncertainty principle but eventually satisfies 
it. Invoking the Peres-Horodecki-Simon's criteria, this corresponds to the 
temporal crossover from an entangled to separable state. The negativity 
and the logarithmic negativity show a monotonic decrease and vanish in- 
dicating the similar crossover. We thus saw that, through the analysis of 
exact dynamics, the effect of environment destroys quantum entanglement 
among Brownian oscillators through the decoherence mechanism. 
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